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ABSTRACT 

Context. The widespread view that cosmological neutrinos, even if massive, are well described since the decoupling redshift z ~ 10 10 down to 

the present epoch by an almost perfectly collisionless fluid of classical point particles is re-examined. 

Aims. In view of the likely sub-eV rest mass of neutrinos, the main effects due to their fermionic nature are studied. 

Methods. By numerical means we calculate the accurate entropy, fugacity and pressure of cosmological neutrinos in the Universe expansion. 
By solving the Schrodinger equation we derive how and how fast semi-degenerate identical free fermions become entangled. 
Results. We find that for sub-eV neutrinos the exchange degeneracy has significantly increased during the relativistic to non-relativistic 
transition epoch at z ~ 10 4 - 10 5 . At all times neutrinos become entangled in less than 10~ 6 s, much faster than any plausible decoherence time. 
The total pressure is increased by quantum effect from 5% at high redshifts to 68% at low redshifts with respect to a collisionless classical fluid. 
Conclusions. The quantum overpressure has no dynamical consequences in the homogeneous regime at high redshifts, but must be significant 
for neutrino clustering during the non-linear structure formation epoch at low redshifts. 
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1. Introduction 

We re-examine the common view in astronomy that cosmolog- 
ical neutrinos, if massive, are well described since their decou- 
pling redshift z ~ 10 10 down to the present epoch by a colli- 
sionless fluid of classical particles (Primack & Gross, 2000). 
Neglected up to now is one of the most fundamental and well 
established principle of modern physics, the fermion/boson 
symmetrisation rules from which Pauli's exclusion principle 
follows (Pauli, 1925; Dirac, 1926; Fermi, 1926). 

When we asked cosmologists how massive neutrinos be- 
have today despite their degeneracy, 1 they would argue that 
they behave as a collisionless fluid of classical particles. In cos- 
mological N-body simulations (e.g., Klypin, et al., 1993) the 
only difference between neutrinos and cold dark matter parti- 
cles is the respective initial velocity dispersion. 

When we asked theoretical and experimental physicists fa- 
miliar with particle physics, quantum and statistical mechanics, 
they would find obvious that degeneracy pressure between neu- 
trinos must be taken into account despite their negligible weak 
interaction cross section, because degeneracy pressure results 
from Pauli's principle and depends only on the exchange sym- 
metry of identical particles. In the limit of negligible interaction 
in the particle Hamiltonian, the fermionic or bosonic behaviour 



1 In the context of cosmological neutrinos, the word "degeneracy" 
is used with different meanings. Here it concerns the particle exchange 
degeneracy, not the mass or the chemical potential degeneracy. 



remains, therefore the particle statistics is an aspect unrelated 
to the strength of the particle-particle interaction. 

In summary, in the cosmologist view neutrinos are well de- 
scribed with a distribution of collisionless classical particles 
following the collisionless Boltzmann equation, while in the 
physicist view neutrinos are quantum particles, i.e., non-local 
entities, and when they are indistinguishable they behave col- 
lectively according to the well established principles of quan- 
tum statistics. In cosmological simulations neutrinos would be 
better represented by a collisional gas with the Fermi-Dirac 
equation of state. 

This distinction has no consequence during the relativis- 
tic homogeneous regime because the initial Fermi-Dirac statis- 
tics is preserved by the Universe expansion (Weinberg, 1972). 
But this is no longer true during the non-linear structure for- 
mation phase, where the different matter components induce 
inhomogeneous perturbations driving neutrinos out of thermal 
equilibrium. Then the degeneracy pressure is different from the 
straight kinetic pressure of collisionless particles, therefore the 
fate of massive neutrinos in galaxy clusters or galaxies differs 
if neutrinos behave as classical or quantum particles. 

In the recent years this problem has taken importance for 
the description of the non-linear phase of structure formation, 
because the likely mass of neutrinos (Tegmark, Vilenkin, & 
Pogosian, 2005) determined in different experiments appear to 
sum up to at least the visible baryon mass fraction. If neutrinos 
behave collectively as a classical collisionless fluid, then their 
interaction with the other matter components is limited to the 
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weak and gravitational interactions, while if neutrinos behave 
as a coherent Fermi gas, then they can develop shocks, increase 
entropy, and lead to non-linear evolution as hard to predict as 
the baryon one during the structure formation phase. 

The paradox raised by the cosmological neutrinos goes in 
fact at the heart of still debated fundamental open problems 
in quantum mechanics, such as the meaning of measurements 
and wave-function reductions, the destruction or built-up of co- 
herence between particles, etc. It concerns also the problem of 
irreversibility in quantum systems, which, as in classical sys- 
tems, are strictly time-reversible in theory, but dissipative in 
practice. 

2. The classical description 

It has been usual in the past decades to assume neutrinos as 
either massless relativistic particles, or as a dark matter candi- 
date able to account for a substantial part of the dark matter, 
which implies a neutrino rest mass well above 10 eV. The clas- 
sical reasoning (e.g., Tremaine & Gunn, 1979; Peebles, 1993) 
goes as follow. Neutrinos are generated above redshift z ~ 10 10 
by elementary particle creation and destruction processes. The 
important point is that at these densities (~ 10 31 cm -3 ) and tem- 
perature (~ 10 10 K) the electro-weak force cross-sections are 
large enough to make the reaction time- scales much shorter 
than the Universe expansion time-scale (see Weinberg, 1972; 
Peebles, 1993; Padmanabhan, 2002). This ensures thermal 
equilibrium and well determined initial thermodynamical con- 
ditions. At redshift below about z ~ 10 10 the creation of 
electron-positron pairs, which releases the last neutrinos and 
anti-neutrinos, drops because the thermal energy of particles 
falls below the 5 1 1 keV electron-positron rest mass energy. The 
neutrino mean free-path for weak interaction diverges and neu- 
trinos are since then considered as propagating freely across 
the Universe, only subject to gravity perturbations. 

The neutrino weak interaction cross section is exceedingly 
small (Peacock, 1999): 

o-~4 - 10 _64 r 2 cm 2 . (1) 

Therefore the widespread view among cosmologists is to con- 
sider neutrinos as collisionless classical particles, i.e., well 
localised mass concentrations moving according to Newton's 
laws. Consequently, neutrinos are supposed to obey accurately 
the collisionless Boltzmann equation. The neutrinos follow an 
incompressible flow in phase space, even during the possibly 
complex non-linear structure formation phase. This leads to the 
well known phase space density constraint (Tremaine & Gunn, 
1979). In this classical description violent relaxation processes 
dilute phase space irregularities into a coarse grained distribu- 
tion, and neutrinos can not increase their degeneracy. 

In the conventional description neutrinos are estimated in- 
sensitive to quantum effects. Peebles (1993, p. 445) argues that 
their de Broglie wavelength 

AlB = — , (2) 

VtlyV 

is a small fraction of cm. The wavelength derived by Peebles, 
0.001cm, is obtained supposing a velocity of v = 300 km s" 1 



and a rather large neutrino mass of m v = 123 eV. Peebles' con- 
cludes: "classical physics gives an excellent approximation to 
the motion of neutrinos, and Liouville's theorem assures us the 
classical orbits will not violate the exclusion principle, because 
the occupation number [...] is conserved". The picture is thus 
that neutrinos are point masses regularly distributed in phase 
space, keeping approximately their relative phase space dis- 
tances sufficiently large to allow neglecting Pauli's exclusion 
principle. A uniformly expanding flow of particles in space is 
a uniformly shearing flow in phase space. Nearby particles in 
phase space do see fluctuating neighbour distances. Therefore 
if any effect depends on phase space relative distances, phase 
space density perturbations can propagate. 

Peebles invokes Pauli's exclusion principle and the de 
Broglie wavelength as criteria for neglecting or not quantum 
effects. Below we will see that relativistic and non-relativistic 
neutrinos are in fact rather degenerate from the beginning, so 
Pauli's principle matters. 

3. The quantum description 

3.1. State of cosmological neutrinos 

Let us check Peebles' assumption with recent estimates of the 
neutrino properties (Dolgov, 2002), keeping in mind that the 
classical world is only an approximation of the quantum world, 
and that the important point is not the ratio of scales (such as 
the de Broglie wavelength vs. the size of cosmic structures) 
but the ratio of the phase space occupation density to Planck's 
constant h cubed. 

In order to not complicate the discussion, we ignore here 
the effect of neutrino oscillations for which the proper collec- 
tive treatment has been only recently worked out (Strack & 
Burrows, 2005). It suffices here to point out that instead of 
discussing the neutrino flavours (the electron, muon and tau 
neutrinos which are superpositions of mass states), we need 
to consider the pure mass states (mi, m2, and m^). The oscil- 
lations introduce coupling terms between the different mass 
states (Strack & Burrows, 2005), and therefore increase the 
coupling between neutrinos species. 

Today's number density of electron neutrinos and anti- 
neutrinos issued from the Big Bang (Dolgov, 2002; Peebles, 
1993, p. 163) of a single family is 3/11 the number of cosmic 
background photons: 

3 * 
ft v ,o = — ft r ~ 1 12 cm . (3) 

Therefore the mean inter-particle distance n~ l J 3 « 0.2 cm is 
comparable to their de Broglie wavelength 

Ad ^ a4cm (iM^r(oT^r- ^ 

So with sub-eV rest mass neutrinos appear rather degenerate. 
Crucial is obviously the value of the assumed velocity. The cho- 
sen v is reasonable for neutrinos bound inside galaxy clusters, 
but it is not obvious whether neutrinos could ever be trapped 
inside clusters. So let us consider the least favourable hypothe- 
sis that identical neutrinos are still homogeneously distributed, 
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Fig. 1. The grey surfaces show 
the entropy per fermion S / kN as 
a function of density and tem- 
perature (log n, log T) for a neu- 
trino mass of m v = 0. 1 eV in 
light grey, and of m v = 10 eV in 
dark grey. The horizontal white 
band marks the constant entropy 
of 4.2018 expected from the hot 
Big Bang fermions. The vertical 
rectangle on the left shows the 
constant density of n = 56 cm -3 
expected for each neutrino state 
at z = 0. The grid lines on the 
entropy surfaces run parallel to 
either constant temperature or 
fugacity. 



with a density of 112 or 56 cm" 3 depending on the still un- 
known Majorana (anti-neutrinos are identical to neutrinos) or 
Dirac (anti-neutrinos are distinct from neutrinos) nature. 

3.2. Entropy, fugacity, and pressure 

The accurate thermodynamical quantities for fermions at all 
regimes can be calculated by evaluating numerically the rela- 
tivistic Fermi-Dirac integrals for particle density n, energy den- 
sity e, and pressure P (Padmanabhan, 2000, p. 216-217) 



n(T,Z) 



e(T,Z) 



P(T, Z) 



4;% r 

h 3 Jo 
h 3 Jo 

^ gs r 

h 3 Jo 



dp 



dp 



Z" 1 exp(e/kT) + 1 

Z" 1 exp(e/kT) + 1 

p 2 1 c 2 p 2 

Z" 1 exp(e/kT) + 1 3 e + rnc 2 



(5) 



(6) 



dp (7) 



where g s is the number of distinct particle states, e = 
^p 2 c 2 + m 2 c 4 - mc 2 is the particle energy, and Z is the fugac- 
ity (the degeneracy parameter), capitalised here to distinguish 
from redshift z. The calculations in this paper have been per- 
formed with the software package Maple, version 10, in which 
the number of significant digits can be chosen. Here 14 digit 
precision has been used. 

The adiabatic expansion of neutrinos since redshift z= 10 10 
means that the entropy per particle S /kN, a dimensionless num- 



ber, has been constant since then. If the initial chemical poten- 
tial is negligible, S/kN amounts to 4.2018 for fermions, and 
3.6016 for bosons (Peacock, 1999, p. 277). The entropy per 
particle can be calculated from Eqs. (5-7) with the relation- 
ship S/kN = (e + P)/nkT, a function of T and Z, which no 
longer depends on g s and h. Calculating the implicit function 
S /kN(T, Z) for a sufficiently wide range of T and Z, we show 
in Fig. 1 its dependence on n and T over relevant cosmolog- 
ical values from high to low redshifts, for two representative 
neutrino masses, 0.1 and lOeV. The grey parts, almost vertical 
in the back, correspond to degenerate states, with n oc T 3/2 on 
the right in the relativistic regime, and n oc T 3 on the left in 
the non-relativistic regime. The almost horizontal parts in the 
front of the S /kN maps correspond to non-degenerate states 
not reached during the Universe expansion. The step passing 
from S/kN = 2.5 to 4 around logT = 2-4 corresponds to 
the transition from the non-degenerate non-relativistic regime 
to the relativistic regimes. The position of the step does depend 
on the particle mass. Clearly, today's cosmological neutrinos 
with sub-eV masses are much colder than the frequently quoted 
1.95 K temperature valid for relativistic neutrinos. 



Observing carefully the constant fugacity curves in Fig. 1, 
one notices that the adiabatic expansion from the relativistic to 
the non-relativistic regime increases the neutrino degeneracy, 
since the constant fugacity curves drop their entropy level. To 
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Fig. 2. Fugacity Z as a function of redshift z for constant comoving 
fermion density n = 56 (1 + z) 3 cm -3 and constant entropy per fermion 
S/kN in the range 4.2018 (lower limit of the bands) and 4.24 (upper 
limit of the bands) for m v = 0.1 eV (lower band), and m v = 10 eV 
(upper band). 



Fig. 3. Ratio of neutrino pressure P to a perfect gas pressure P c = nkT, 
as a function of redshift z for constant comoving fermion density n = 
56 (1 + z) 3 cm -3 and constant entropy per fermion S/kN in the range 
4.2018 (lower limit of the bands) and 4.24 (upper limit of the bands) 
for m v = 0.1 eV (lower band), and m v - 10 eV (upper band). 



make this point quantitative, we solve numerically for T and Z 
the pair of equations 

^-(T,Z) = 4.2018, n(T,Z) = 56 cm" 3 , (8) 
kN 

with a bisection search technique because the non-linear 
equation solver of Maple for floating-point number solutions 
f solve based on a Newton-Raphson method fails in this case. 
The value of n is a low value for today's neutrino density 
(choosing ft = 1 12 cm -3 would reinforce our conclusions). We 
find that the temperature scales as the inverse of the mass, while 
Z at z - is practically independent of the neutrino mass. We 
obtain for the present neutrino temperature, 

r y0 = 1.614 • 1(T 4 (—) 1 K, and, Z = 16.50. (9) 
VleV/ 

So between the creation of neutrinos at z ~ 10 10 and now, 
somehow fugacity has increased from 1 to 16.5. 

To show how fugacity behaves as function of redshift at 
constant entropy per particle, we ask Maple to solve numeri- 
cally by interpolation the above quantities in Eqs. (5-7). The 
result is shown in Fig. 2. Fugacity increases from the early 
Universe to the present by a factor 16.5: massive neutrinos 
are substantially more degenerate now than at their decou- 
pling epoch. This is independent of their mass if they are non- 
relativistic. When calculating the Z(z) relationship for a range 
of entropies, we see that if for any reason entropy would in- 
crease instead of being strictly constant, fugacity would in- 
crease even more. 

Associated to fermion degeneracy is the Fermi pressure 
contribution. To quantify this important factor for dynamics, 



again we ask Maple to evaluate numerically the integrals in 
Eqs. (5-7), and to find by interpolation the pressure at constant 
specific entropy curve. In Fig. 3 the ratio of the full pressure to 
the perfect gas pressure P c = nkT is displayed. The additional 
quantum pressure with respect to a collisionless gas amounts to 
5% at high redshifts, independent of the neutrino mass. At the 
transition from the relativistic to the non-relativistic regime, the 
redshift of which does depend on the mass, the quantum pres- 
sure contribution increases to 68%, and becomes subsequently 
independent on the mass in the non-relativistic regime. Again 
any increase in specific entropy would increase pressure. 

Obviously this additional pressure is a substantial effect at 
low redshifts determining the clustering of neutrinos during 
the structure formation phase. The immediate prediction is that 
neutrinos should tend to accumulate less in bound structures 
than truly collisionless particles. However, dissipative effects in 
non-linearities such as shocks may eventually falsify this pre- 
diction. 

3.3. Asymptotic estimates 

These previous results are exact down to the used 14-digit nu- 
merical precision. Let us now also derive intermediate states 
with asymptotic approximations, in part to estimate the errors 
made with this more traditional approach. 

During the relativistic phase of expansion neutrinos cool 
proportionally to the photon temperature (Peebles, 1993): 

/ 4 \ 1/3 / 4 \ 1/3 

Tv ' r = \nj Ty = \n) To(l+z) ' (10) 
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where T Q is today's CMB temperature of 2.73 K. When the neu- 
trino thermal energy equals its rest mass energy, neutrinos be- 
come non-relativistic. Then 



m vr c = 3kT v 



4\" 3 

it) : 



3£r (i+z). 



Then the critical temperature is 
Tvr = Hh^* 3970 



3k 



and the critical redshift 
l^-19«(^). 



(ID 



(12) 



(13) 



Below redshift z c temperature drops for adiabatic non- 
relativistic matter approximately as 

T " = T " (r+zf * 98 • 10-4 <> + *" (tjv)"' k ■ (>4) 

Thus the thermal speed of cosmological neutrinos amounts to 



' J_\ 1/3 3kT 
11/ m v c 



(l+z)« 151 



kms _1 .(15) 



(see Ringwald & Wong, 2004, for a similar derivation). With 
respect to the exact result Eq.(9), T^ n is overestimated by a 
factor 6.07, and v v by a factor 2.46. 

Owing to number conservation, the neutrino density n v 
varies in proportion of the volume 

n v = n v , (l +z) 3 . (16) 
The average phase space volume per neutrino is estimated as 



1 4ttvJ ^ U47T(kT 

V Yly 3 11^0 \ C 



(17) 



which is independent of the neutrino mass and redshift, as long 
as neutrinos are non-relativistic. Compared with the elementary 
quantum volume h 3 , we have, 



Vy _ 284 

h 3 n y? o 



5.1 



(n v ,o \~ 
V 56 / 



(18) 



Neutrinos appear therefore rather close to be strongly de- 
generate, consistently with the previous accurate calculations. 
Therefore quantum effects must be considered with more atten- 
tion. 

4. Quantum Effects 

4.1. The paradox of collisionless identical particles in 
the quantum regime 

In view of the previous results, we are led to a paradox. On 
one side from cosmological conditions we find that neutrinos 
are fairly densely packed in phase space, so Pauli's exclusion 
principle should be applicable, while on the other hand in the 
usual view among cosmologists neutrinos should ignore each 
other since their particle-particle electro-weak coupling is very 
low. How can particles without interaction term interact? 



Often in quantum mechanics one must use the idealisations 
of isolated, independent systems, although wave-functions for- 
mally extend over the whole space. This is justified (e.g., 
Cohen-Tannoudji et al., 2000, chap. XIV, D, p. 1384) when the 
system wave-function is itself well separated from the outer 
world wave-functions. This is possible because particles are 
emitted not as plane waves extending uniformly over the whole 
space, but as wave-packets since any emission process takes a 
positive time-interval. Therefore the individual wave-functions 
consist of well defined square integrable functions with com- 
pact support, as required by quantum theory. 

So if we picture individual neutrinos as localised moving 
wave-packets, they follow a classical free particle trajectory 
as long as the wave-packets are well separated. But when the 
wave-packet cores of identical particles overlap, then quantum 
interference must occur, blurring the possibility to distinguish 
the particles, because they must simultaneously follow Pauli's 
rule of anti-symmetrisation for fermions. Pauli's symmetrisa- 
tion rules for fermions and bosons are independent of the par- 
ticular particle kind and strength of the interaction potential. 

All the interactions between elementary particle, nuclei, 
atoms, molecules, and even larger pieces of matter like super- 
fluids, follow the Pauli symmetrisation rules: whenever par- 
ticle wave-functions overlap, quantum mechanics is consis- 
tent with experiments at the condition of applying the wave- 
function symmetry postulate. Somehow the individual packets 
of identical particles which could be considered as indepen- 
dent when far apart, must form a symmetric or antisymmet- 
ric combined wave-function when they overlap. The speed of 
this (anti-)symmetrisation process has nothing to do with the 
electro-weak interaction, but just the quantum mechanical na- 
ture of identical particles. 

Before discussing this paradox further, let us illustrate 
graphically with exact solutions of the Schrodinger equation 
how the wave-packets of two identical free fermions behave. 

4.2. Interaction between two identical free fermions 

It is not difficult to solve exactly the time-dependent 
Schrodinger equation, 



ifi -¥(x, t) 
at 



HVix, t), 



(19) 



for the Hamiltonian H of two identical free fermions and illus- 
trate graphically the behaviour of the most physical object rep- 
resenting a particle, the amplitude of the wave-function core. 
We illustrate below what happens to two spinless fermions in 
3D space without explicit interaction potential. As initial con- 
dition we adopt Gaussian wave-packets because they allow to 
calculate all the results in closed form. 

As derived in standard quantum mechanics, the full so- 
lution of the single free particle Schrodinger equation in 3D 
Cartesian space reads (with units normalised conveniently for 
our cosmological neutrino context to [/] = 1 cm, [m] = 0.1 eV, 
and [t] = 1.7 • 10" 7 s, which allows to set h/m = 1): 



1 



(Imf) 3 ' 2 



x ¥(y,0)d 3 y. 



(20) 
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Fig. 4. Interaction of two free fermions in the classical regime starting initially in opposite directions and horizontally (k x = 6, k y = 0, k z = 0). 
Shown is the normalised squared norm -x, t)\ 2 in the plane (x, y, 0). 



Thus finding how an arbitrary initial wave-function ^(x, 0) 
evolves is just a matter of calculating a triple convolution in- 
tegral at any wished time t. For particular initial ^(x, 0) such 
as Gaussian packets the integral can be completely expressed 
in closed form, which we manipulate below with Maple. 

When describing two independent distinguishable parti- 
cles, the full wave-function takes the form 

¥(xi,x 2 ,0 = ¥i(*i, f)¥ 2 (*2,f), (21) 

for which the time evolution is just the product of the single 
particle solution Eq.(20). 

But this is only a very particular form of all the wave- 
functions x ¥(x\,X2,t) that the Schrodinger equation for two 
particles can describe. Functions that cannot be separated as 
a product of single particle wave-function describe entangled 
states. In particular two identical particles must follow either a 
symmetric or anti- symmetric wave-function upon the exchange 
of the particles. These read 

x ¥ ± (x u x 2 j) = 

-L[w l (x u t) x ¥ 2 (x 2 , t) ± ^(xz, f)V 2 (x u t)] , (22) 

where the + sign applies to symmetric functions and the - sign 
to antisymmetric functions. The time evolution of such an en- 
tangled pair of wave-functions can also be calculated exactly in 



closed form by forming the product and sum of single particle 
solutions. The explicit formula is too long to show here, but can 
be easily manipulated with Maple or similar software tools. 

Although the wave-function of two particles is 6- 
dimensional and complex, when the particle centres of mass 
and velocity are set at the origin we can focus the attention on 
the function ^-(x, -x, t) since any measurement of a particle 
at location x will constrain the position of the other particle 
at position -x, since we assumed to know the centre of mass. 

-x, t)\ 2 describes the probability of finding one particle 
at position x and the other one at -x. 

Below we show graphically the normalised value 
-x, t)\ 2 in the plane z = for the initial condition in the 
centre of mass and velocity frame of two fermions represented 
initially by a Gaussian wave-packet of unit half- width: 

¥i(x, 0) = exp(iifc • x) exp [-(x - x ) 2 /2] . (23) 

The initial position of particle 1 is Xi(0) = {-5, -2, 0}, and the 
initial momentum is k\ = {k x , 0, 0}. The initial position of par- 
ticle 2 is symmetric about the origin in order to put the centre 
of mass and velocity at the origin. The momentum is therefore 
all in the component k x , which is not directed toward the centre 
of mass in order to illustrate a generic interaction with non-zero 
impact parameter. 

Here we show three different representative values of the 
momentum k x illustrating the classical (k x = 6), semi-classical 
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Fig. 5. As in Fig. (4), but in the semi-classical regime (k x - 2) 

(k x = 2), and quantum regimes (k x = 2/3). This should suffice 
for our purpose of illustrating the exchange interaction of free 
identical particles. 

In Fig. 4 we see the two wave-packets individually spread- 
ing slowly due to their high momentum k x = 6. Despite that, for 
a while the packets strongly overlap, subsequently they recover 
their localised wave-functions. In such a case the individual 
packets behave as expected for classical particles: they move in 
straight line and preserve their identity. Although Pauli's prin- 
ciple applies, it has no practical consequence asymptotically. 
These two particles can be approximated as distinct again at 
later times. 

In Fig. 5 we see that at lower momentum, k x = 2, two 
wave-packets interfere much more strongly. The antisymmetric 
wave-function always vanishes at the mass centre. The size of 
this exclusion region around the origin grows at lower veloc- 
ities. The interference pattern presents elongated wave crests 
which are inclined with respect to the original particle momen- 
tum directions. A plane wave seems to be building-up in a di- 
rection slanted with respect to the particles directions. Not only 
the positions are strongly mixed in a single structure, but the 
momentum becomes also mixed. After a while the individual 
wave-packets representing the classical localised particles suc- 
ceed to partly reform, but a remaining coma of entangled states 
remains around the origin. 



At still lower relative momentum, k x = 2/3, quantum ef- 
fects become dominant. In Fig. 6 the two particles initially 
distinct rapidly loose their identity and form a single anti- 
symmetrical wave-packet which spreads around the origin with 
waves again reflecting an intermediate direction of propagation 
with respect to the initial particle momentum. Clearly in this 
case the probability of finding a neutrino somewhere else than 
along the initial straight trajectories is strongly enhanced by the 
quantum interferences. The classical approximation of straight 
trajectories is no longer valid, and the initially distinct particles 
are entangled forever. In the quantum regime a clear dispersive 
behaviour occurs with marked asymmetric past-future asymp- 
totic states. 



The same experiments have been repeated for bosons. The 
only notable difference is that the wave-functions are always 
positive at the origin, and even take there maximum values in 
the quantum case: this corresponds to the Bose-Einstein con- 
densation. In contrast, the fermionic wave-functions always 
vanish exactly at the origin, which means that a volume of 
space around the origin is hardly available to the particles. 
Fermi pressure is a consequence of confining the particle en- 
ergy in a reduced volume. 
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Fig. 6. As in Fig. (4), but in the quantum regime (k x = 2/3) 



4.3. Interaction between N identical free fermions 

The problem of describing statistically N identical bosons or 
fermions with or without interactions has been solved very 
early in the history of quantum mechanics (Uhlenbeck & 
Gropper, 1932), but is much more difficult to illustrate graphi- 
cally than the 2 particle case. What is clear is that in the limit 
of vanishing particle interaction, the bosonic or fermionic sta- 
tistical behaviour tends toward classical Maxwell-Boltzmann 
statistic only at high temperature. At sufficiently low temper- 
ature the statistical behaviour tends toward the Bose-Einstein 
and Fermi-Dirac statistics respectively. 

An interesting aspect developed for example by Huang 
(1964, chap. 10), is that the first quantum correction to a classi- 
cal perfect gas, which is appropriate for our cosmological neu- 
trino context, caused purely by the bosonic or fermionic na- 
ture of the particles is mathematically equivalent to an particle- 
particle interaction potential, with the peculiarity to depend 
also on temperature T, 



<p(Xi - Xj) = -JcT log 



1 ± exp -2n 



(24) 



Fig. (7) shows that the pseudo potential is repulsive for 
fermions and attractive for bosons. The potential becomes im- 
portant ((p ~ kT) at a distance slightly below the de Broglie 
thermal wavelength /L dB = h/ ^InmkT . 



—\ — n 1 1 1 1 1 1 1 r 



T 1 1 1 1 1 - 




Fig. 7. Effective particle-particle interaction potential for a perfect 
gas of fermions (top curve, repulsive potential), or bosons (bottom 
curve, attractive potential) in first order of quantum correction, as a 
function of distance r in unit of the de Broglie thermal wavelength 
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This is consistent with the rule that quantum effects become 
important when thermalised particles, are less distant in the 
classical description than A^b • The potential permits to derive 
an effective interaction cross-section, and from there the time 
required for particle to become entangled. We will refer to this 
time as the coherence time-scale. If the particles were classical 
we would call this a collision time-scale, that would determine 
a relaxation time. 

If we adopt a classical description of neutrinos moving in 
such a potential, clearly a strong deflection of trajectory occurs 
whenever (p/kT ~ 1, which means also that \x t - Xj\ ~ A^b, 
so the interaction cross-section reads cr ~ 7rA^ B . Thus the first 
order quantum correction to a classical description of semi- 
degenerate neutrinos introduces an interaction potential pro- 
ducing "entanglement" with a characteristic time- scale given 
by 

1 ^ n%' 2 y/kr 
TiyCTVv n v h 2 

. 13 . (25) 

using the asymptotic relations for T, n, and v given in Sect. 3.3. 
This yields a very short time-scale «c 10" 6 s for any reason- 
able present day and past cosmological neutrinos. Although 
this derivation is not consistent since quantum particles do not 
follow classical trajectories, at least r characterises the entan- 
glement time-scale between neighbour neutrinos. 

4.4. Decoherence time 

In the decoherence theory (Zeh, 1970, 2005; Zurek, 2002, 
2003), classical physics emerges from the much larger sensi- 
tivity of the entangled states in the density matrix x ¥(x) x ¥*(x'), 
x ± x' , to the interferences coming from the outer world 
than the diagonal states x ¥(x) x ¥*(x). Particles become local- 
ized classical objects because the uncorrelated perturbations 
from the outside world prevent the wave-function to spread. 
The decoherence time To for a quantum system considered 
over a distance Ax to loose quantum correlations x ¥(x) x ¥*(x f ), 
Ax < \x - x*\, is related to the relaxation time r R induced by the 
outer world, and to A$& ' 




Thus the decoherence time is inversely proportional to the 
square of the considered region. Obviously, to quantify deco- 
herence we must identify the fastest process that perturb or "re- 
lax" neutrinos. 

Let us make a best effort to find the shortest possible tr that 
cosmological neutrinos might be subject to. If we discard weak 
interaction processes (cf. Eq. (1)), we only have gravitational 
deflection available. 

Suppose we have a density n of bound gravitational objects 
of mass M. The gravitational cross-section for neutrinos at ve- 
locity v v is cr = nG 2 M 2 /v v ~, and the timescale for deflecting 
significantly neutrinos by these objects is: 



nervy nG 2 nM 2 



Thus at constant mass density nM more massive objects deflect 
faster than light objects. Galaxies seem the most promising de- 
flectors, because stars are mostly bound in galaxies and amount 
to only a fraction of their mass, while galaxy clusters contain a 
small fraction of the galaxies. 

Plugging in plausible numbers for galaxies (n « 
10" 74 cm" 3 , M « 10 45 g), and v y « 10 8 cms _1 , we obtain 
tr ~ 10 22 s. Therefore for td to be less than the entanglement 
time < 10" 6 s, we need to consider regions larger than 

Ax = A dB ^H L ~ 10 13 cm, (28) 
V r D 

in order to begin to use classical physics. In such a region of AU 
dimension we have ~ 10 40 neutrinos that behave as a fermionic 
ensemble. Since the relaxation time ~ 10 22 s is much larger 
than the neutrino age ~ 10 17 s, most neutrinos that would be lo- 
calized by hypothetical detectors would still be entangled with 
neutrinos at cosmic distances apart. The ratio of these times, 
10" 5 , gives the fraction of neutrinos that have been perturbed 
by the gravitational attraction of a galaxy. After 10 5 entangle- 
ment times of < 10" 6 s, each neutrino has a large probability 
to be also entangled with perturbed neutrinos, thus propagating 
the coherence loss. 

5. Discussion 

5.1. The status of quantum physics 

The quantum physics rules are so strange that it took a long 
time for the physicist community to make a standard view 
about its principal aspects. Today several points are not com- 
pletely resolved, such as the meaning of the measurement pro- 
cess, or the role of time in a relativistic quantum theory of sev- 
eral particles (Penrose, 2005). Quantum gravity is still in the 
work. 

However in the last decade many new experiments have 
confirmed the strange rules of quantum mechanics. The 
Einstein-Podolsky-Rosen (EPR) paradox has been confirmed, 
two entangled particles can fly apart over long distances and 
preserve their mutual correlation alive even when one of the 
particle is subject to measurement. In recent experiments (Tittel 
et al., 1998) the photon entanglement has been preserved over 
distances of the order of 10 km throughout a complicated net- 
work of telecommunication optical fibres. 

5.2. Hanbury Brown and Twiss (HBT) effect and the 
quantum decoherence 

At the root of quantum optics, there is a fine example of ex- 
change symmetry interaction operating on otherwise interac- 
tionless particles: photons. This effect was first investigated by 
Hanbury Brown and Twiss (Hanbury Brown & Twiss, 1954) 
with the now famous intensity interferometer (Hanbury Brown, 
1974). The effect could also be performed at the laboratory 
scale (Morgan & Mandel, 1966) and was an important step in 
the development of quantum optics. In the Hanbury Brown and 
Twiss experiment, independent incoherent photons emitted all 
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over the surface of a star (D ~ 10 9 m) at distinct times are col- 
lected on Earth by two distinct telescopes separated by macro- 
scopic distances ( d ~ 10 2 m). The collected distinct photons ar- 
riving in the two detectors within a photon coherence time are 
verified to be correlated in time. Despite being particles with- 
out explicit coupling term and initially incoherent, photons suf- 
ficiently close in phase space (they are very close in momentum 
space) succeed to bunch, i.e., to exchange momentum. This is 
an example of a partial Bose-Einstein condensation. 

The fermion effect corresponding to the HBT experi- 
ment has been confirmed recently for electrons in laboratories 
(Henny, et al., 1999; Kiesel et al., 2002). Such experiments are 
complicated by the fact that the repulsive electro-magnetic in- 
teraction term of electrons must be cancelled by either an atom 
lattice, or a sufficiently large distance between electrons. In 
these experiments, electrons overlapping their wave-functions 
are found to be anti-correlated, which is at the root of the de- 
generacy pressure in fermionic gases. Better suited particles are 
neutrons, for which the HBT effect has also been verified re- 
cently (Iannuzzi, et al., 2005). 

The important point of these verified phenomena is to show 
that identical particles interact also through their exchange 
symmetry, and that the speed of this kind of interaction, not ex- 
plicit in the Hamiltonian, is in practice determined by the speed 
at which the core of the wave-functions overlap significantly. 

If we take the strong overlap of the particle wave-functions 
as the principal requirement for particles to be said interact- 
ing, and the weak overlap of many uncorrected wave-functions 
as the regime during which the decoherence due to the rest of 
the world acts, we can understand why our common experi- 
ence of the world is not a fully entangled structure (Penrose, 
2005), but the classical world that we know, where spatially 
distinct regions and classical devices can be assumed to ex- 
ist. Correlations between particles are constantly lost by deco- 
herence and built by wave-function overlaps. The macroscopic 
result for multiparticle systems is equivalent to a relaxation to- 
ward a thermal state (Gemmer et al., 2004). 

Therefore if we apply these considerations to cosmologi- 
cal neutrinos, these must overlap their wave-function very fre- 
quently (<sc 10" 6 s), and each "interaction" produces a partial 
entanglement of the wave-functions. Quickly neutrinos get en- 
tangled to a growing number of other neutrinos which together 
form a spreading entangled wave-function. Any perturbation 
or "observation" of a given neutrino entangled with others con- 
cerns then all the other neutrinos. Everything happens as if the 
detected neutrinos had exchanged momentum with the others, 
like classical particles in usual interactions do. 

6. Conclusions 

In view of the solidity of quantum mechanics and the likely 
sub-eV mass of neutrinos, we arrive to the conclusion that 
quantum effects play an dominant role for cosmological neutri- 
nos. Despite being almost interactionless, these neutrinos form 
a non-localisable ensemble of entangled particles subject to 
Pauli's principle. The weakness of their interactions with the 
outer world is precisely the key factor that prevents them to 
decohere faster than they become entangled. The only way to 



justify a classical description of neutrinos following the col- 
lisionless Boltzmann's equation would be to find a perturba- 
tive process that would decohere all of them in much less than 
10" 6 s. But then this perturbative process should presumably be 
included in the neutrino dynamics. 

With the most accurate description of particles, the wave- 
function, we have illustrated explicitly how a pair of initially 
localised neutrinos becomes entangled fast due only to their 
antisymmetry; quickly the localised wave-packets become a 
single fully non-localised object. In practice this means that 
more neutrinos act rapidly as an ensemble that must be treated 
with the standard quantum statistics tools. Since, except for an- 
tisymmetry, no interaction potential is significant, the appro- 
priate description at cosmological scales is simple, the one of a 
semi-degenerate perfect fermion gas. 

For structure formation, the Fermi pressure contribution 
has little importance as long as the Universe is homogeneous, 
since only pressure gradients play a dynamical role. At red- 
shifts z » 10 5 the Fermi over-pressure amounts to 5% of the 
kinetic pressure. But at redshifts z 10 4 in the non-relativistic 
regime, during the non-linear structure formation epoch, the 
over-pressure increases to 68%. Then a Fermi gas of massive 
neutrinos is substantially different from a collisionless gas. Not 
only the velocity dispersion remains then isotropic in a tightly 
coupled gas, but since collisional fluids are able to shock due 
to the non-linear convective term in the fluid equations, the 
evolution of the combination of neutrino gas, dark matter, and 
baryons in the today Universe may turn out to be complex. 
Shocks typically increase entropy, degeneracy, and therefore 
pressure, as seen in Sect. 3. It would be interesting to deter- 
mine how far neutrinos could violate the phase space density 
conservation by dissipative processes, especially when mixed 
with baryons. 

The extension of these results to other dark matter parti- 
cle candidates resulting from the Big Bang should be straight- 
forward. Relativistic fermions, degenerate or not, can not con- 
tribute significantly to matter clustering. Usual heavy cold dark 
matter fermions, such as neutralinos, are assumed much less 
dense and much heavier than neutrinos, in which case phase 
space degeneracy may remain small and decoherence fast, al- 
lowing the use of the collisionless fluid approximation. 
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